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Abstract 

We present a comprehensive theoretical description of quantum well exciton- 
polaritons imbedded in a planar semiconductor microcavity. The exact non- 
local dielectric response of the quantum well exciton is treated in detail. The 
4-spinor structure of the hole subband in the quantum well is considered, 
including the pronounced band mixing effect. The scheme is self-contained 
and can be used to treat different semiclassical aspects of the microcavity 
properties. As an example, we analyze the "selection" rules for the exciton- 
cavity mode coupling for different excitons. 
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I. INTRODUCTION 



Since the first experiment on quantum well (QW) excitonic polaritons in the semicon- 
ductor microcavity (SMC)S extensive experimental and theoretical studiesE have been con- 
ducted on excitonic polaritons in QW's and multiple QW's embedded in an SMC, and even 
bulk excitonic polaritons in an SMC. The SMC exciton-polariton was first described by so 
called linear dielectric modelEl which was originally proposed for atoms in a microcavity. 
Since then, most analysis on the microcavity problem has been more or less based upon the 
assumed analogy between excitons and atoms. The non-locality of the dielectric response of 
the exciton in a QW and the intricacy of real excitons due to the complex structure of the 
valence band has not been treated seriously. 

In fact, excitons differ from an atomic excitation in that they are a sort of excitation 
from electrons in the valence band confined by the QW's barriers. The excitonic polariton is 
essentially a quantum many-body effect and can be understood to be a result of the medium's 
polarization. In the calculation of the dielectric response to light propagation in a QW, we 
have translational invariance only in the plane perpendicular to the growth direction, so 
the nonlocality of the dielectric response along the growth direction should be carefully 
treated. This makes the problem of a QW exciton embedded in an SMC non-triviallii. 
In particular, the barrier confinement effect in association with the 4-component spinor- 
like hole wave function should appear in the nonlocal dielectric response. The confinement 
induces hybridization between heavy and light hole subbands, which sometimeJHi,0 plays a 
nontrivial role at the top region (F-point) of the valence band and has interesting observable 
consequences for the optical transitions of the excitons in a QwiHi. This makes the non- 
local dielectric response of the QW more complex. To our knowledge, the non-locality of QW 
exciton-polariron has been consideredS for the simplified exciton model, while the real QW 
exciton model including non-locality and subband mixing effects has not yet been considered 
seriously, especially in a microcavity environment. On the other hand, although the effect 
of such hybridization induced exciton is not very strong, the exciton-photon coupling is 
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greatly enhanced in the microcavity. Thus it is interesting to investigate the nature of such 
hybridization induced QW excitonic polaritons in the SMC environment. 

Based upon the above considerations, in this paper we provide a self-consistent semiclas- 
sical description of the excitonic polaritons in a QW embedded in an SMC, in which the 
hole subband hybridization, the nonlocality induced by the QW barrier confinement and the 
boundary conditions for the SMC in connection with distributed Bragg reflectors(DBR's) are 
all consistently taken into account. The light field is treated semiclassically but the polariza- 
tion of the medium in the QW is treated in the context of quantum many-body theory with 
the electromagnetic wave-electron-hole interaction, so that the description of the effective 
coupling between light and excitons does not need any phenomenological input "coupling 
constant" . The only input parameters are the width of the exciton levels. The calculated 
Rabi splitting of conventional excitonic polaritons could be considered as an improvement 
to the so-called dielectric model in which the photon-exciton coupling constant is estimated 
using the electric field of a cavity mode and and the oscillator strength of an exciton. We 
will show also the complexity in introducing such an effective oscillator strength due to the 
coupling of the motion along or perpendicular to the growth directions for the hybridization 
induced polaritons. As an interesting effect of the barrier confinement induced electron and 
hole envelope wave functions for a system with a pair of symmetric DBR's, we apply parity 
symmetry analysis along the growth direction to the DBR's and the SMC confined QW, 
and obtain a kind of selection rule for excitonic polaritons in the well, which is also quite 
useful for finding the hole subband hybridization induced polaritons. As an application of 
our scheme, we have calculated the reflection spectrum for several different excitons and for 
different incident angles, and have obtained some interesting results, which will be published 
elsewhere. Here we only discuss the formal aspects of our approach. 

Since this topic is in an interdisciplinary field between condensed matter physics and 
quantum optics, we present our discussion in a self-contained way for readers in both fields. 
This paper is organized as follows. In Sec. II., we formulate the description of the exci- 
tonic polaritons in a QW embedded in an SMC in terms of a propagating electromagnetic 
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field, while the SMC is further confined by a pair of DBR's. The basic spirit of its classical 
electro dynamic aspect mainly follows ref.il, but the crucial difference is in the physics con- 
tained in the formal expression for the dielectric response. Reference^ is valid chiefly for the 
intra-band transition while our paper is devoted to the interband transition, particularly the 
hybridization effect of the hole band. Sec. Ill is devoted to a quantum many-body descrip- 
tion for the medium polarization in a QW, in which the hybridization for the heavy hole 
subband and the light hole subband is emphasized. Although the complex valence band in a 
QW and the related exciton problem has been widely treated in the literatures, a complete 
treatment that starts from the quantum many-body theory for the QW exciton, and includes 
non-local response function and complete symmetry analysis for the 4-component spinor like 
wavefunction, is still absent. We shall present such a complete and detailed treatment in this 
section. In Sec. III-A, we give an expression for the medium polarization propagator with 
both positive and negative frequency parts, which is applicable to generic non-translational 
invariant systems. In Sec. III-B, by summing over the Kramer's degenerate states, we show 
the detailed expression for the nonlocal conductivity tensor which is proved to be diagonal 
and can be factorized into a sum of bi-products. This makes its non-locality effect explicit. In 
Sec. IV we present the final complete set of self-consistent equations as well as the boundary 
conditions which are put into a discrete version for the convenience of practical calculation. 
As an application, we discuss the "selection rules" for the exciton-cavity coupling. Sec.V 
gives some concluding remarks. 

II. DESCRIPTION OF A QUANTUM WELL IN A CAVITY CONFINED BY 
DISTRIBUTED BRAGG REFLECTORS. 

The SMC under consideration consists of a QW of thickness A embedded in a thin layer 
semiconductor material which is sandwiched between a pair of distributed Bragg reflectors. 
This semiconductor layer serves as the barrier for the electrons and holes in the QW, while 
it forms an optical cavity confined by the DBR's. The whole structure is schematically 
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illustrated in Fig. [l|. Although we only treat the case for one QW in the cavity, generalization 
to more complex cases can be easily done. The thickness of the cavity is with £c(^) as its 
medium's dielectric constant, where uj is the light radiation frequency. We assume that the 
QW medium has the same (background) dielectric constants as those of the cavity medium. 
The left (right) DBR is constructed from A^^. {^r) pairs of two alternating layers. One of 
the two layers has a length Li and dielectric constant ei{uj) while the other L2 and £2(^^)5 
respectively. We choose the the growth direction with z = being the center of 

the system. Thus, the QW is located at (— f-, f ) and the SMC occupies a region from — ^ 
to with Lc > A. Consequently, —kf {kf) is the boundary between the left (right) DBR 
and the SMC. Moreover, we rotate the coordinate axis around the 2;-axis in such a way to 
make the incident light propagate in the x — z plane. Then the wave vector has the form 
q= {q^ = gil, 0, Qz = q±) with Qy being always equal to zero. As a result, for the p-polarized 
wave, the electric field has only x— and z— components, and Ez, respectively, while for 
the s-polarized wave the electric field has only a y-component Ey. For a propagating wave 
with fixed frequency and qj] = q\\ex (where e^, fi = x,y, z is the unit vector along the /i-axis), 
its electric field has the form 

E{f, t) = E{ij, qii , ^)e'('?ii— *) = E{uj, rje"'"*. (2.1) 
The Maxwell equation for the electric field can be written as 

2 A 

^ (V' + - VV] ■ Eiu, T) = -^J>, r-) (2.2) 



in which "i" is the medium index. Taking into consideration Eg.( |2.1|) , Eq. ( p.2|) can be 
simplified to 



Li{uJ,q\\) ■ E{uj,q\\,z) = —j{uj,qi\z) (2.3) 



with 



= - N (2-4) 

For simplicity, from now on we discuss mainly the p-polarized wave propagation. It is 
straightforward to convert the discussions for a p-polarized wave into those for an s-polarized 
wave. We will do this when necessary. Outside the QW medium, the polarization current 



density j{uj, 2) = in Eq. (p73|) . We then have plane wave solutions with 
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For the left DBR, denoting the x-components of the incident and reflected amplitudes of the 
electric field as Aq and Bq, which are defined as approaching to and leaving from the right 
boundary of the left DBR. By applying the method of ref.0, we obtain the corresponding 
electric field amplitude at the right boundary of the left DBR as Ai and Bi, which are also 



defined by approaching the boundary from the left as 

nL rpL 
11 5 -^12 



^0 \ / Til 1 \ / ^1 



-DO/ \-'2l5 -'22/ V-Dl 

with 

\ ^21 ' ^22 / V ^21 ' '^22 / 

where is the transfer matrix of the left BDR, is the tranfer matrix of a pair of layers 
with different dielectric constants, as a unit cell of the left DBR, in the DBR and is the 
number of the pairs in the DBR. Then for the right DBR, we introduce further ^3 and S3, 
the forward (along the z-direction) and backward (in the negative z direction) amplitudes 
of the electric field, respectively. They are defined by approaching from the right to the left 
boundary of the right DBR. Meanwhile, A^ and B^ are the corresponding amplitudes at the 
right boundary of the right DBR which are also defined by approaching the boundary from 
the right. We can similarly define the transfer matrix by solving the homogeneous Maxwell 
equation^! to obtain 

(2.8) 



(2.6) 



13 \ / J 11 , 12 \ / ^4 

bJ [ri, Ti) U4, 



with 



' rpL rpR \ / ,R J.R \ Nji 

11 5 12 \ / ''11 1 '-12 > 



rpR rpR. / I ,_R ,R 

-^21 5 -^22 / \ ''21 1 '•22 . 



(2.9) 

The detailed expressions for the matrix elements tf^- and tf^ with i,j = 1,2 are shown in 
Appendix A. 

In the QW region, there is a sort of additional polarization effect of the medium other 
than that described by its background dielectric constant. This is because the virtual 
electron-hole pairs creation and annihilation processes will renormalize the light propagation 
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in the QW region. In particular, the Coulomb interaction between the virtual electron-hole 
pairs may play a crucial role in a certain frequency range. In this section, we just introduce 
formally a conductivity tensor a{uj, qf, z, z') to describe this effect, such that 

= J dz'a{u,qfz,z') ■ E{uj,q\\,z') (2.10) 
which will be further investigated in the next section. Therefore, in the SMC region (in- 
cluding the QW), the Maxwell equation (|2.3| ) and (|2.4| ) can be converted into an integral 
equation as 

E{u, q\\- z) = E(^\u;, gi,; z) - ^ J%dz' J% dz"G{u, gji; z, z') ■ 

a{uj,q\\-z',z")-E{LO,q\\,z") (2.11) 

in which 

i(^)(a;, q\\]z) = iAj^^'^' + B^e-'^^-^')e^ - ^iAj''^'^^ - B^e-''^^'^')e, (2.12) 

Q±,c 

is the homogeneous solution of Eq. (|2.3|) , while G{lu, ^j; z, z') is the Green's functionil for the 
equation 

L{uj, gJi) ■ G{uj, qf z, z') = 5{z - z'). (2.13) 
We can easily solve Eq.( p.l3|) up to a homogeneous solution of itself with its explicit expres- 
sion shown in Appendix A. We notice that, in all equations, from Eq.( 2.1l| )to Eq. (|2.13|) , z 
and z' are confined to the region 

- — < z < — (2.14) 
2 - - 2 ^ ' 



Moreover, the integral equation (|2.11|) exhibits the nonlocality not only provided by the 



Green's functions, but also induced by the medium polarization effect which is described 
by the conductivity tensor. To solve this integral equation, we emphasize that it should be 
solved with proper boundary conditions (BC's). This is because neither the homogeneous 
solution E^^\uj, q^^, z) nor the Green's functions are solved with respect to the correct BC's 
for the SMC. Therefore, following classical electrodynamics, at 2; = ±^ we should have the 
EC 

Ai + Bi = A^e-''^^'-^ + Sce*«^.-^ 

rA/2 r / rA/2 i ii (-< Y, , ^, . _ -f-c 
f.2 l^v'=x,z l^i>"=x,z J-A/2 ^xy\^^q\\i 2 ' 



X au',u"{uj, q\\;z', z' 



)Ey„{uj,q\\;z") 



(2.15) 



- Bi) = ^{A^e~"^^'-^ - S,e^<?^.-^) 

+ 'fzV^ Y.v'=x,z J2u"=x,z /-A/2 dz"Gxy (w, % 




) 



(2.16) 



Si/'=x,2 Siy"=x,z /_A/2 /-A/2 G^x,;/' (^1 "2^! ^ 

X (T;,',,,"(u;,^|; 2;',2;")E^//(a;,gJ|; z") 



) 



(2.17) 




(2.18) 



for fixed cj and since Aq is the input while B4 is usually taken to be zero. Here we 
have altogether eight independent constants 50,^1,-61,^0-80^35-83 and A^, and a two- 
component electric field function Ex{uJ,qi\, z) and Ez{uj,q\\, z) as unknown variables (func- 
tions), in which E^ and Ez are defined only in the SMC (including the QW). We stress 
that they are mutually coupled to each other through Eqs. ( p^.6| ), (^.8|),( prT3| )-( [^.18| ) and the 
integral equation Eq.( |^.llD (referring also to Eqs. ( |2.12| ) and (|XT|)-([7^). Eqs. (P^ ) and 
( ^78|) are 2x2 matrix equations, the above equations are exactly eight mutually coupled 
algebraic equations and a two-component integral equation, which solves self-consistently 
the whole DBR-SMC(QW)-DBR system as long as we have the detailed expression for 
the conductivity tensor. Such a description^ constitutes our mathematical framework for 
studying semiclassically the excitonic polaritons in the QW embedded in a SMC . The pair 
of DBR's is coherently correlated with the SMC and plays the role of mode selection. Thus, 
the existing modes in the SMC not only dynamically couple to the semiconductor medium 
in the QW, but also coherently couple to the DBR's. This is the physical meaning of such 
a description. 
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III. POLARIZATION OF THE SEMICONDUCTOR MEDIUM IN A QUANTUM 

WELL 



In this section we concentrate on the polarization of the semiconductor medium in a 
quantum well. As mentioned in the introduction, the loss of displacement invariance in the 
growth direction of the QW and whole microstructure makes the dielectric function nonlo- 
cal. The non-locality of the medium polarization is closely related to the wave functions, 
and these functions are closely related to the band structure of the semiconductors. Thus, 
first of all, we should give a brief description of the band structure of the semiconductor 
microstructures. It is well known that the electronic states in a microstructure such as a 
QW can be satisfactorily described by the envelope function approximation. However, in 
many publications, this approximation is often an oversimplified version, i.e. the electron 
and holes are described as free particles with simple effective masses. Such a simplification 
usually works well in most cases, especially for conduction band electrons. However, this 
is not true for valence bands, i.e. for holes. The top of the valence bands of the semi- 
conductors with Td symmetry belongs to the Fg irreducible representation, which is 4-fold 
degenerate. In bulk material, the valence states apart from the F point can be divided into 
heavy and light holes, both of which exhibit strong anisotropy and non-parabolicity. For the 
microstructure case such as in a QW, the situation is more complex. Generally speaking, 
any valence wavefunction can only be described by a four-component spinor, the basis of this 
spinor i.e. the Bloch cell periodic function (with spin) satisfies certain space transformation 
relations, and more importantly, for different components the envelope functions have dif- 
ferent space rotational transformation properties. This invalidates the simple free particle 
picture. Although the problem of hole subband structure and related exciton's as well as the 
optical transitions have been treated by several authors the complete analysis based 
on the exact non-local dielectric function and four-spinor nature of the wave function is still 
absent. We will devote the following section to this problem. 

For the usual III-V compound semiconductors, at the conduction band bottom Tq, the 
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Bloch cell periodic function u^^\f) is a spin doublet with index a = |, — |, while at the top 



of the valence band Fg the cell function u^^) (r ) can be regarded as a four component spinor 
with a' = |,|,— |,— |, the HH subband corresponding to a' = ±|, and the LH subband to 
a' = ±|. The plus and minus signs refer to time reversal (Kramer's)degenerate states. In 
the effective mass approximation, the non-interacting Hamiltonian for the envelope functions 
written in the first quantization representation has the form 

1 9 \2 ^ 1 /h d Y 



;/3,/3' 



n. 



a,f3 



Eg - ^ + Ve{r\le,Ze) 



(3.1) 
(3.2) 



zrh 



1 

2m 



(Pi, Q, 



R*, 
V 0, 



0, 



R, 
0, 

A, 



R*, Q, 



\ 
R 

-Q 

Pi J 



(3.3) 



where 



A = (7i+72)(f4) +(7i-272)(f^ 
F2 = (7i-72)(f4)' + (7i + 272)( 



n a 

i dzh 



= -«V3o.(fi)((f5f,)-<f4)) 

2\ 



h a 

i dVh 



In Eqs. 



i-(|3.4|), the conduction electron coordinate is giving by fe 



1,6) 



(3.4) 
Zp) with 



f^l_e = XpCx + yeGy, and the valence hole coordinate by = (rj|,/i, Zh) with = XhC^ + yh(iy; 
Eg is the gap which separates the conduction band and the valence band, /i the chemical po- 
tential, m the physical electron mass, and and ml are effective masses of the conduction 
electrons corresponding to motion perpendicular or parallel to the z-axis; and are the 
confinement potentials which form the barriers for the QW. For the DBR-SMC(QW)-DBR 
system discussed in this paper, Ve and Vh are actually rj|^e and rjjjj independent. Eq( p.3|) 
with ( |3.4|) is the well known Luttinger Hamiltonian^! which is written down in a form with 
a specific choice of the coordinates. In these two equations, Q and R describe the hybridiza- 
tion between the HH and LH subbands, while 71, 72 and 73 are band structure parameters. 
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We notice that the envelope function description adopted in this paper is vahd around the F 
point. We denote ff^'^ir) envelope function which can be solved from the eigenvalue 

equation for the Hamiltonian equations ( |3.2| ) and (|3.3| ) with appropriate boundary condi- 
tions in connection with the confinement barriers. Its corresponding energy eigenvalue is 
denoted by ef'^\ In this notation, b = c refers to the conduction band with A = ±| corre- 
sponding to a spin doublet; h = v refers to the valence band with indices A = ±|, ±| which 
describe valence band branches resulting from diagonalization of the hybridized Hamiltonian 
Eq. (|3.3|) . Notice that each of the eigenfunctions (envelope functions) corresponding to these 
branches is a four component spinor. We again name the envelope functions with A = ±| 
as the HH branch and that with A = ±^ as the LH branch, according to the properties that 
the A-spinor has the dominant component as a = A. Moreover, s is the quantum number 
depending on the confinement potential. It is discrete along the direction of confinement 
but continuous in the extended directions. 

A. Semiconductor Medium Polarization Induced by Virtual Electron-Hole Pairs 

It is known that the conductivity tensor is connected to the medium's polarization tensor 

by 

afo;; r, r ') = — Tfrfcu; r, r ) (3.5) 

UJ 

Intuitively, during light propagation, the dominant contribution to the medium's polariza- 
tion should be the virtual excitations of electron-hole pairs. To clarify certain conceptual 
problems we shall discuss the ideal situation: the intrinsic state of the QW at low tem- 
perature, i.e., the valence band is almost fully filled while the conduction band is almost 
completely empty. The virtual pairs thus consists of conduction electrons and valence holes. 
In such a case, the Coulomb interaction should induce a series of exciton states distributed in 
the semiconductor gap. These are bound states formed by conduction electrons and valence 
holes. These virtual bounded e — h pairs should also contribute to the medium's polariza- 
tion even at zero temperature for an intrinsic QW. Based upon such an understanding, the 
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polarization can be calculated straightforwardly as outlined in Appendix B. We notice that 
there are no requirements for the spatial translational invariance, so that it can be applied 
to any sort of QW. This also has the advantage that when applied to the planar QW, the 
in-plane center of mass momentum for the excitons can be explicitly treated, which is a 
non-trivial property for the excitonic polaritons. Now, the derived expression has the form 
nriuj;f,r') = (^)^EaEa' J2pEi3'{v,a'\p\c,a){c, /3\p\v,a) ■ 



'/'L/ (^.^)V' '{r,r) 



_^„(A,f;V.s';n) f) (3-6) 

in which 

where V'ia'*''*' '"'*(^ei "Oi) is the two body exciton-like wave function for the ri-th eigenstate 
which consists of a conduction electron with quantum numbers A, s in (rg, a) representation 
and a valence hole with quantum numbers A' and s' in the {fh, a') representation. It satisfies 

= i?(^'^^^'^Vit'^''^'"\r-;,f,), (3.8) 
where the eigenvalue E^^'^'-^' '^''^ is exactly the energy term appearing in Eq.( |3.6| ). In Eq. (|3.6| ), 
we have also introduced the dipole matrix element 

{c,a\p\v,a') = 1 / d'r^{ut\f)W:){f) - V<W(f) ■ u^:}{f)). (3.9) 
where f2 is the volume of the crystal cell and the integration is also confined within the 
cell. We emphasize here that, for the inter-band processes, the dipole matrix element is 
carried by the Bloch cell functions but not by the envelope functions. In Eq.( |3.6| ) only the 
two-body wave functions with its arguments Te = = r contribute to the spectral weight 
of polarization propagator, r is the point where the electron meets the hole . Moreover, the 
first term in the bracket of Eq. ( p.6[ ) is its positive frequency part while the second is the 
negative frequency part. Both parts are necessary for a boson-like propagator. We notice 
that the r and r ' change into each other in the two corresponding spectral weights, for 
which we understand that the positive and negative frequency parts correspond to a sort of 
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forward and backward propagation, respectively. Only when a certain kind of symmetry is 
present for the system under investigation will the two spectral weight functions be equal 
to each other. This is one of the main topics for the next sub-section. 

We notice that the extended eigenstates of Eg. (|3.8D can also contribute to the spectral 
function. Since they are off-resonance from our interested frequencies, their contribution 
has been included in the background dielectric constants. 

B. Summation over the Degenerate States Based on Space-Time Inversion 

Symmetries 

After performing a summation over the space-time reversal transformation connected 
degenerate states, we show in this subsection that not only the positive and negative fre- 
quency parts in the polarization propagator can be combined into one term, but also the 
non-diagonal elements for the polarization tensor can diminish. This will make the expres- 
sion neat and further simplifies the calculation a great deal. 

We apply the formal discussion in the previous sub-section to the system we interested 
in - a planar QW embedded in a DBR-SMC-DBR system (see Sec. II). The general subband 
indices and quantum numbers A, s; A', s' and n are now specified into the following symbols: 
gjj; Hex, I] n,jc] n', and jh- Here gj| is the 2D center of mass momentum for the exciton, which 
is a good quantum number describing the translational invariance in the X-Y plane for the 
system; n^x and / are the quantum numbers for the eigenstate of the exciton describing the 
in-plane relative motion for the virtual electron and hole pairs; Uex is the major quantum 
number while / is the angular quantum number; n {n') is the index of the discrete states 
due to the confinement in the z-direction of the conduction electrons (valence holes); jc = 
±^ are the subband indices describing the spin doublet for the conduction electrons while 
j/i = ±|, ±1 are the subband indices for the valence holes representing the total angular 
momentum induced by spin-orbit coupling. 

To solve (fe, fh) from Eq.(|3.8|) we need the complete set of eigenfunctions 
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of the non-interacting Hamiltonian described by Eqs ( p.2|) and ( p.3|) for the conduction 
electrons and valence holes, respectively. These are 

^e^^1i--1i-V^(r)(ze)5,.,„ (3.10) 

for the conduction electrons and 

(3-11) 



for the valence holes with subband hybridization, where is the extension for the QW in the 
X-Y plane. In solving Eqs. (|3^ ) and (|3.3|) , there is an energy eigenvalue associated with 



Lp[^\ze) for the z-direction electronic motion and another eigenvalue e^'^j^ = = 0) 

associated with (p^^)^j^.^,{k\\^h-, Zh) which could also be interpreted as describing the ^-direction 
hole motion. Actually, for a hole in the valence band, its motion along the ^-direction 
described by the eigen wave function (Pn']jh;a' coupled to its motion in the X-Y plane 
characterized by k\\^h- So (Pn\jh;a' depends not only on Zh but also on k\\^h whilst ^n'ljh ^^^^ 
has the k^^h dependence. Such a subtlety is due to the hybridization term Q and R in 
Eq. (pp and (gj) 

We notice further that, in the QW region, vl^" \'^\\,h, z) is constant so that we chooce it 
to be zero. The 2;-axis parities of the matrix elements of the Luttinger Hamiltonian Eq. ( p. 3D 
change their signs alternately along each row and each column. This interesting symmetry 
property makes the hole eigenfunction for Eq. (|3.3| ) have an amazing parity symmetry so that 
the first and third components (a' = |, — |) have the same parity symmetry along the z-axis 
which is opposite to that of the other two components (a' = Moreover, since the 

Luttinger parameters 72 and 73 are usually rather close to each other, we may assume that 
they are equal, i. e. 72 = 73. We further introduce a polar coordinate {kw^h^'&h) for k\\^h in 
the X-Y plane, k\\^h = k\\^h{cos'dh€'x + sm{}hey). Following Eq. (|3.4| ), we have Q ~ e"*'^'* and 
R ~ e^^*'''". As a result, the corresponding eigenfunction in Eq. ( p.llD acquires the following 
functional dependence 

^u'LA^U^ = ^'^''""'^''^uLAh^^h) (3-12) 

in which -dh is separated from the product k\\^hZh is a single dimensionless 
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argument. 

We then expand the two-body exciton-hke wave function in terms of the non-interacting 
electron and hole wave functions, Eqs. ( |3.10D and ( |3.11D for a fixed set of good quantum 



numbers q\\ ; Hex, I', n, jc, n' ,jh, and obtain 

(g||;n,,,/;nj-n',i;,) ^ ^ e^(^||,e-n|,e+fc||,„-r||,^) . 

&'\he, kh) ■ ¥'i^H^e)5,.,«^S,,^„,(^i|,., Z,) (3.13) 

The expansion coefficient /^"^"'"""H^IU) can be understood as the 2D excitonic wave 
function for which the center of mass degree of freedom has been separated. If we ignore 
the hybridization terms Q and R in the hole Hamiltonian Eq. ( |3.3|) , and further ignore the 
Ze and Zh dependence for the Coulomb interaction in Eq. (|37^ ) because the thickness of the 
QW is rather small, we can show easily that f^J^i'^''\k\\^e-,k\\^h) would be precisely the 2D 
hydrogen atom solution of Eq. (|3l8| ) constrained by ^| = fcy^e + Since now we have the 
hybridization, the physical meaning of /^"^" ''"'''(^11,6) is no longer so transparent, and 
the coefficient should be determined by substituting itself into Eq. ( p.8| ). Actually we need 
only the value of the wave function ( p.l3| ) at re = Hi = r. Thus, we have 

^11 ,e ^11 ,h 

fS'''\ke, k'^)^l^\^)^nLAku^ ^) (3-14) 

Moreover, the GaAs system has an energy gap ~ 1.5 eV. The wavelength of interest has an 
order of magnitude ~ lO^A. On the other hand, the effective Bohr radius for the 2D-exciton 
in GaAs has the order of magnitude ~ lO^A. Therefore, k\\^e and k\\^h "will span a region 
in momentum space in which the 2D excitonic wave function /^"'" '-"'H^lhe) has non- 
negligible contributions, much bigger than that of ^| . We may reasonably set q\\ to be ~ in 
the double summation for A^y e and k\\^h-, i-e., we may have approximately k\\^^h — ~k\\,e = ^|| 
in Eq.( p.l4| ) within the double summation. Then 

= 7pe^^''-"'''^..,«/(£f./S"''^(fc||)v^i^)(^ (3-15) 



Intuitively, the 2D exciton-like wave function should have the expression 
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/£?^^(A^||) = /S"'^(^ll)^''' (3-16) 
where {k\\,i)) are the polar coordinates of A;||. Then, by further utihzing Eq. |3.12| and 

Jq'^ d'dexp i{l + jh — a')'d = 27r5;+j,^ we obtain 

= lo^ ^/i3'"'^^(^ll)v'l^H^)¥'i'k;«'(^lF)- (3.17) 

We substitute expression p.l7| ) into Eg. (|3^) with its indices changed from the general form 



into the present, and obtain 



/^ _ e2 1 



roo fc||rffc|| roo ^||'^^f| f Jh) Z L \ f (";"•' Jh) ^ N 
Jo 2tt jo 2n Jn^x,l \'^\\)Jn^:,,l \'^\\) 

"P^a^ {z)V^n^ (Z') ¥'Sj^;a'=i;,+« (^11 z) ■ ( A;|| z') 

where 

E^'^'im) = E, + eS;^"'^'') + e'^^^^^ + 6^ + e^?,^ (3.19) 
depends on q^^ = only. In Eq. (|3.19|) , e^"'"''-^'') is the kinetic energy for the center of mass 
of the 2D exciton and ^n'^iJ'/-"'^ < is its binding energy. The e^*^^ and e'l^) are associated 
with the confinement potential for the conduction electrons and valence holes which has 
been introduced in connection with Eqs. (|3.10|) and ( p.ll|) . As mentioned above, due to the 



hybridization of the valence holes, it is difficult in practice to reduce E^^^{-''^\q\\) into such a 
neat form as Eq. ( p.l9| ) with conventional interpretations. The function f^J^i^''\k\\) used in 



Eq. (|3.18|) is essentially the eigenfunction of the eigenvalue problem of Eq. (|3.8| ) which can 
be solved either approximately or numerically. In practice, a variation solution was used. 

Besides the dynamical calculations, we have various symmetry properties for the system 
which are exact. Due to the space-time reversal symmetry, i?^"^" (gy) and /l"^" (^11 ) 
are degenerate with respect to ±/, jc = ±|, jh = ±| and jh = ±|. We have also the time 
reversal symmetry property for the hole wave function 

^^:LA^) = ^iU-,-A^)- (3.20) 

Furthermore, all the f^J^y^^\k\\) ■, Vn\^) and V^i'^j^;Q,'(fc||-2) are real functions. On the other 
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hand, it is known that there is only one independent matrix element for {v,a'\p^\c,a), i.e. 
p = {X\px\s) = {Y\py\s) = {Z\pz\s) where {X\, {Y\ and {Z\ are the orbital part of the 
Bloch cell functions for the hole band with / = 1 and \s) is the orbital part of the Bloch cell 
functions for the conduction band. By making use of the Kramer's degeneracy properties 
and the Clebsch- Gordon coefficients for the matrix elements {v, a'|p^|c, a), we can sum over 
the plus and minus values of l,jz and jh but keep the absolute values of /, jc, jh and the 
corresponding energy eigenvalue fixed. After a lengthy calculation, the conductivity tensor 
associated with Eqs. (|3.5|) and (p.l8|) can be derived in a diagonal form as 



(j^^(u;;f,r") = -^^e"'"'^"'!! a^^{uj; qf z, z') (3.21) 

911 

with 

(^tiui^', q\\', Z, z') = ^^Sfj,^u En En'.ifc J2n^,,l E/^ 

n) ■ 0SAH-)e^'''^v££rH^'), (3.22) 

in which 

<^E;/:^(^') = f (3.23) 

2 p(n;n',jh)/ \ 

p{n;n',jh) _ ^^n^x,l /„ 

In Eqs. (|]22D,(|]23D and ( ^24]) , jh takes the value of | and \ for HH and LH bands. 



respectively, while / is equal to and 1, //i = if / = 0, and//i should be summed over ±1 

if / = 1. This "new" quantum number 1^ results from the summation over the degenerate 

states due to Kramer's and space inversion symmetries. The term rj'^^h\iM comes from the 

Clebsch- Gordon coefficients for the dipole matrix elements, with 

1 



„{3/2A0) ^ „(1/2,1,1) 
//i Ijj. 







j^a/2,0,0) ^ ^(^3/2,1,1) = ^a/2,1 -1) ^ 

77(3/2,1,1) =0. (3.25) 
We stress that if we ignore the hybridization between the HH and LH subbands, then only 
/ = exciton states contribute to the spectral weight function. Meanwhile, ^n'',jh;a'=jS^\\'^'^ 
becomes Lp^^}j^.^,^j^ (z) which is independent of k^\. Therefore, the hole subband hybridization 
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induced excitonic polaritons are characterized by the angular quantum number I ^ 0. We 
stress further that, in Eqs. ( p.22| ) and ( p.24| ), we replace u + ir] hy u + ^7!"^" ''''''' which 



7n"l" the only phenomenological parameters introduced in this approach for describing 
the width of the exciton level. 



IV. DISCRETE REPRESENTATION OF ELECTRIC FIELD IN AN SMC 

In the last section, we obtained the polarization of a medium induced by the Coulomb 
interacting virtual electron-hole pairs in an intrinsic semiconductor QW at zero temperature. 
We notice that the z, z' dependence for the conductivity tensor Eq. ( p. 221) is separated into a 
bi-product of two (/)-functions. Thus, the kernel of the integral in equation ( p.ll| ) is separable. 
Following Ref.ii , we can easily transform the integral equation (|2.11|) for the electric field 
in an SMC into an algebraic equation which could further simplify the calculations. In 
particular, as Pn^^i ''^'^\uj, q\\) has a resonant pole, only very few discrete components for the 
electric field are needed, and possibly only one may be sufficient. 

Introducing 



= j_^/^€:sr:\^)E{u^A-^-)^ (4-1) 

we may express the electric field in terms of E^^^^{j'^\uj , q\\) as 

E(.;,gl,2) = EW(u;,gl,^) + ^x 

En,n',i^ T^n^^Uu ^^^'^S''^" V, % , Z, z') 

with 

^{j„iM = S^^nl^'^'^M (4.3) 

We emphasize that the F-functions introduced by Eq.( |4.1| ) depend on the electric field de- 
fined only in the region of [— |-, but the corresponding electric field Eq.([4.2|) is meaningful 
for the whole region of [— -y]. Utilizing the above two equations, we can express Eq.( p.llD 



as well as the boundary condition Eqs. ( |2.15D -( p.l8D in terms of F-functions. We perform 
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the operation J^j^^^ dz(j)^2:%ii''\^) both sides of Eq. and introduce f{uj, 



A/2 
-A/2 



'^{n;n',jn;n";n"',j'^). 



dz(f), 



/-A/2 /-A/2 



A/2 



-A/2 



A/2 



and 



(4.4) 



(4.5) 



from which we obtain 



B{ri;n' ( , _ f{n;n' / , \ ^ 

ti,{n;n',j,; n";n'",j'^) {n";n"',j'J 

■V^^''''^-F^^;fiu^.%)- (4.6) 
The boundary condition Eqs. ( p.l8|) can also be reformulated in terms of F^^i.f^^'^ as 

^J- En.n'jh ^ne^,l,lh /-^A/2 ^z' P^^^'^^i '^''\uj, q\\)G x,f,{uJ , q\\\ 

^ m^c2 ^n;n',jh ^nex,l,lh /-A/2 ^zT^^^ i '^''\uj, q\\)Gx,ij. 
A3 + B3 = Ace^«^-^=/2 ^ Q^^~iq^,,Lj2 



^ 



(4.7) 

(4.8) 
(4.9) 



V Sn.n'jft T,ne^,l;lh /-^A/2 ^zP^^'^^i '-"'^ (t^, ^| ) G^^,.,^ (w, g|| ; ^, 2^') 



^(^3 - S3) = -i£-(A,e^''^'=^=/2 ^ ^^g-igx,cLc/2) ^ 



(4.10) 



^2 T ,p(n;n',i,J 



J2n;n',jh J2nex,l,lh /-A/2 '^^ -^71ea;,Z 



(4.11) 



Then, in terms of the discrete representation for the electric field F^^^i.i^\u , q\\) , the com- 
plete description of the DBR-SMC(QW)-DBR system is now given by Eqs. (^-( |4.11|) and 
(|;6D, referring also to Eqs.(g;|), (PD , (p:T2|) , ([X3D , (^) , (|3)and ( p^ . Note that in order 
to express /i"!^" ;z'^''(i^, Q\\) in Eq. (|4.6| ) in terms of Ac and Be we need to introduce two more 
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terms w ^^^^'ifi^ (o;, ^| ) and w' ^^^'^'^^^ (u;, gj, ) 

= j'^'^^/^^~''^-''&:\^)E^'\^ (4-12) 
according to Eq |2.12| . Up to now, it is clear that we do not need the effective photon-exciton 
couphng constant any more in the approach we have adopted. 

We have now estabhshed a complete set of equations which can give us the real space 
electric field distribution within SMC with non-local dielectric response of QW excitons em- 
bedded in it, and thus enable us to obtain various semiclassical optical properties. Compared 
with previous theoretical work, our approach treats both the electro-magnetic field and the 
electron excitation exactly and does not depend on any artificial parameters. Moreover , 
certain effects not considered in previous studies such as violation of simple selection rules 
and the strange angle dependence of the reflection spectrum, can be easily calculated using 
this approach. We will report these results elsewhere. As an application, we analyze some 
interesting aspects of the exciton "selection rule" for a symmetric microcavity. 

It is well known from ordinary exciton models that only the S exciton can couple to 
the electromagnetic field, but due to the four component nature of the hole subband wave 
function and the existence of band mixing, the p exciton can also couple to the photon field 
and thus form polaritons in a microcavity. We call this a hybridization-induced exciton- 
polariton, and analyze the relevant parity propeties along the z axis. 

It can be seen that Eq. ( |4.2| ) in conjunction with expressions (|3.23|) and ( |4.1| ) shows 
explicitly how the non-locality for the conductivity behaves. Each term of the spectral 
weight function has the form of a product of two 0-functions with one of the two depending 
only on z and the other on z' . Therefore, one of them is convoluted with the electric field and 
plays a role as a "weight" function while the other is convoluted with the Green's function 
as a "source" function, which strongly influences the effective coupling between light and the 
exciton. This makes the role of the parity symmetry along thr 2;-axis explicitly important. 
Here we restrict ourself to the case of normal incidence and symmetric DBR pairs only. 
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As we have discussed in Sec. II-B, v^.[f^(^) and ip^^'^j^^.^i{k\\z) can be classified by the parity 
symmetry (along z-axis). In fact, the parity symmetry for the contributive component of 
the hole spinor wave functions is controlled by its index a' = jh + Ih- Moreover, for an 
SMC confined by a pair of symmetric DBR's, it is known that E^'^^lu , q^, z) as a function 
of z is also an eigenstate of the parity symmetry at the resonant condition Lc = |A, with 
n = 1,2, ■ • • . By considering Eqs. (|CT| ), {\K^ , and (g]^, we can verify that the 

electric field E(lo, z) will have the same even-oddness as that of E^'^\lo, q\\; z). Following 
Eqs. ( ^^) , (p.4|) and (|3.23| ), the physics contributed by the QW depends strongly on the 
integral 

= l_^^^dzcp(:\z)^^:}^^.^^,^^^^^^ (4.13) 

and we can now find interesting consequences for the polaritons. If we set ip!;^\z) for the 
conduction electron to be always in the n = 1 state, it will not contribute to the parity 
symmetry since it is an even function of z. For a symmetric SMC, E^^\ij,qi\; z) is even for 
Lc = ^ and odd for Lc = A. Hence for an SMC of Lc = |-, the electric field is even. Then, 
following Eq. (|3.17|) , only polaritons with quantum number I = Ih = and the index a' = jh 
can survive in the SMC for both jh = \ and jh = |. This is because ip^^)^i^j^.^i ^j^{k\\z) is an 
even function of z. The hybridization induced polaritons with quantum number / = 1 cannot 
be observed since t/:'^'^Lij,^;a'=j^±i(^||^) are odd functions of z which will make ^^^.^^ = 

0. On the contrary, if we have an SMC of = A, the electric field Ex{uj , q^f, z) becomes 
an odd function of z, if (p^^)^-^^j^.^,^j^{k\\z) is still an even function, the polaritons with 
quantum number n = n' = 1, I = will be forbidden for the A-SMC which will make 
Ii^'^'''''\u, q\\; k\\) = 0. In this case, only the hybridization induced polariton can survive 
in the SMC since (^^'^Lij^;Q'=j^±i^(^|j^) is odd. What we learned from the above discussion 
is that the even-oddness of the electric field depends on the cavity resonance condition, 
while the parity of the components of the electron or hole wave function depends only on 
the QW. Then we can have an interplay of the even-odd symmetry for the integrand of 
Eq.( [4.1| ). As a result the parity symmetry along the z-axis will provide a sort of selection 
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rule for the forbidden polaritons in a symmetric DBR pair confined resonant SMC. This 
could be helpful for finding the hybridization induced polaritons. Our calculation shows 
the existence of HH subband dominated polaritons which have the quantum numbers n=l, 
n'—2, jh—\i'nex = 2, I — 1 with splitting value ~ 0.2 mev. The details and other interesting 
results will be published elsewhere. 

V. CONCLUDING REMARKS 

In conclusion, we have presented a self-consistent semiclassical approach for exciton- 
polaritons in a QW embedded in an SMC. In this approach, the effect of the complex 
valence band structure and the non-locality of the dielectric response of the exciton in the 
QW are carefully considered. For IHH excitons and normal incidence, our approach gives 
the same results as those obtained previously. For complex cases such as high index excitons 
and oblique incidence we expect that it could predict some new phenomena. For example, 
we have shown that a 2P exciton can also couple to a photon mode and form a polariton. 
Moreover our analysis gives a "selection rule" for the formation of exciton-polaritons in a 
symmetric SMC, which is essentially an interplay among the angular quantum numbers of 
excitons, the electron-hole subbands indices and the resonance conditions of the SMC. 
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FIGURES 



FIG. 1. Schematics of the symmetric planar semiconductor microcavity 



FIG. 2. The summed Coulomb interacting electron- hole bubble diagram series 



APPENDIX A: THE TRANSFER MATRIX ELEMENTS AND GREEN'S 

FUNCTIONS. 

The transfer matrix elements can be calculated straightforwardly a si 



and 



''11 ~ ^ 



cos(g^,L,)-i(|g + ||^)sing^,,L, 



''12 ~ ^ 



''22 



''21 



-ig±,iLi 



i f £2g±,i siq±^\ ■ T 

- 77 ■ — sm o I 2-tvs 

i ( ^2g±,i 

.2 VeiiJ±,2 



£ig±,2 
£29^,1 



-) sing_L,2-^^2 



cosg,,L,-|gf^ + |f^)sing,,L, 



''11 ~ ''ii ' 



fR _ y-i 
''22 ~ ''22' 



'-12 



'-21 ; 



''21 ~ ''121 

Moreover, the Green's functions in Eq. (|2.13|) can also be easily solved aS 



G^^^{u}, gy; z, z') = G,^^{uj, qf z, z') = ^jf (u.) ^^ 9n{z - 2;')e^'?x,ck-^'l 



Gz,z{^^,q\\]z,z') 



2 

_£! §(z - z') ^gjgx,c|2-^'l 



Gy^y{uj,qi\;z,z') 



i ^ipi r\z-z'\ 



Gxy Gyx Gyz Gzy 
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APPENDIX B: OUTLINE OF THE DERIVATION OF THE POLARIZATION 

PROPAGATOR 



We outline here the main steps for deriving Eq. by a summation over the Coulomb 
interaction ladder diagrams accommodated in a bubble diagram consisting of a conduction 
electron and a valence hole, as shown in Fig. ^ 

We introduce the second quantized wave operator and current operator: 

m = EEEE¥'&'nr)«(f)(r)a^^) (Bi) 

b X s a 

j{f) = -^{ij^{f)Vip^{f) - Vip^{f)'4j{r)) (B2) 
2imc 

in which the second quantized electron annihilation operator for a state with band 

index c, A and quantum number s. Note that we dropped the index k of the cell-periodic 
functions of the Bloch function, in corresponding with the envelope function approximation 
. Using the linear response theory approximation and the Mastubara representation, the 
polarization part of the electron system 7rr(iu;; r, r ') can be expressed as 

nriuj; f, r') = - / dre*-"^Tr{/)T,j(r, T)j(r o)}|..„^.+., (B3) 
Jo 

where p is the thermal density matrix, Tr the trace operation, and the chronological 
operation along the imaginary time axis. Expression (p^) has the advantage that it can be 
calculated systematically by applying the diagrammatical technique. We notice first that 
each pair of electron and hole lines in the upper series of diagrams in Fig. ^ contributes a 
term 

where s' describes the quantum numbers charge conjugated to s'. Each pair of electron and 
hole lines in the lower series of diagrams in Fig. § with the directions of the arrows reversed, 
contributes a term 

in which ^^'^''^^ = e^^'^^ + Eg — fi, ^^'"''^'^ = e^^''^^ + fi and n(^) is the Fermi distribution function 
for the conduction electrons and valence holes. Moreover, a dotted Coulomb line contributes 
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a factor 

-^i-^^' = E [Mr^:^^in^'::'i^)-T^^^^ (B6) 

where Sc = Sc{uj) is the background dielectric constant for the medium. Finally, each ver- 
tex at the end point of the bubble contributes a dipole matrix element {c,a\p\v,a') or 
{v, a'\p\c, a) . All the above expressions follow straightforwardly from the quantum many- 
body text-book with an additional consideration shown in (|Bl| ), i.e. our basis wave function 
is not the usual simple plane wave, but the envelope function associated with a more "mi- 
croscopic" cell periodic function u. 

It should also be noticed that, for an intrinsic semiconductor in the low excitation limit, 
the conduction band is almost completely empty and the valence band almost completely 
full, so the two Fermi distribution functions in Eqs. (p^) and (p5|) can be taken as zero. Fur- 
thermore, this guarantees that only the ladder diagrams give any contributionll. Through 
summation over the Coulomb interacting ladder diagrams accommodated in a bubble dia- 
gram which consit of a conduction electron and a valence hole (shown in Fig. we can 
obtain the following expression for the polarization tensor, after a direct but lengthy calcu- 
lation, 

7r{iuJn; r, f') = (^)^ Ea,a' Ea,. Ex' ,s' {vta\r)ip*}^^^ \f) {v , a'\p\c, a) 



-vtc^'\ncp^::'j\f){c,a\p\v,a') 

1 

. -.(c,A) ^(v, 



(B7) 



^si'i'^^n, r') = Ea",a"' '/'sfa'^ '/'i'.a"' (^') (c, a"\p'\v, a'") 

, f 2\ A.A":A"',A' 
+ ) l^\"^s" l^X"',s"' 'f^s,s":s"',s' ■ 



1 



T:„ :,„ '{iuJn\ r') 



ll^n — t, II —t,-iii 



(B8) 

/ 2\ Y^ A",A;A',A"' 

~ \^ ) l^X",s" l^X"',s"' ^s",s;s',s"' 
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We then introduce two auxiliary functions 



(BIO) 



and 



(Bll) 



We can express T[^^!'\iuJn,r') and f r') in terms of the two auxihary functions by 

utihzing the orthogonal properties of ^pi'^aK'^ ^^"^ f^s''a'\^^ i'^ ^ similar to an inverse 
generalized Fourier transformation in the Hilbert space. We have further the completeness 
relation for ^'f^^^r) and f^^'^a'^^')- With the aid of these procedures, the polarization tensor 
can now be derived as 



',r') = {^yj:a,a'{{v,a'\p\c,a)P, 



- (c, a\p\v, a')Pa,a'{iuJn; r, r; f')}. (B12) 
Next, we apply the operator iiOnSa,a"Sa',a"' — [Ho]a,a';a",a"' to the Pa"^a"'{i^n',f,f'',r') from 
the left, apply the operator iuJnSa",aSa"',a' + [Ho]a",a"';a,a' to Pa",a"'{i^n] f,f']r') from the 
right, and obtain 



\^iiL!mSa,a"Sa',a"' — (^Hq 



. ' i dr' ' i df' 



\f-r'\Ja,a';a",a"'} 



Pa" ,a"' {iuJm; r,r'; r') = {c,a\p\v,a')6{r - f')S{r' - f') 

J2a" Ea'" Pa",a"'{'i'^^m; f,f'] r')\iUTn5a",a^a"',a' + (-^^0 VCd.~l 1 d 



i dr ^''^ ^ i 'df' 



\f—r'\ 



-) } = {v,a'\p\c,a)6{r ~f')6{r' ~f'). 



(B13) 



We may then solve the auxiliary functions from Eqs. ( [B13D and ( [B13| ) by the Green's 
function method to get 

Pa,a'ii'^m;r,r;f') = J2Ga,a';p,p'ii'^m]r,f;f'f'){c,P\p\v,P') (B14) 

with 



X,s X',s' n 



p(A,s;A',s') 



(B15) 
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as well as 

Pa,a'{i(^m;r,f;r") = '^{v,P'\p\c,f3)Gi3,i3'-,a,a'{iuJrn-,r,r-,r,r) (B16) 

/3,/3' 

with 

/ (A,s;A',s';n) / -»/\ ; *(A,s;A',s';n) /- ~/\ 

Gp,f3';a",a"'{tuJm; r , r ; r , r ) = }^ }^ }^ ■ ^(x,s;X' ,s') 

A,s X',s' n liOm + rjn ' ' 

in which "0^^/ ' {f,r ■') and are exactly those introduced in Eg. p.^). By sub- 

stituting Eqs. (|BT^)-(|BT^) into Eq.QBT^), we obtain Eq.(^ as desired. 
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Fig.l of Xu et. al.,"Semiclassical. 
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Fig. 2 of Xu et. al.,"Semiclassical..." 





c:e conduction electron 
v:h valence hole 
Coulomb interaction 



